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Abstract 

We argue that in the measured px domain of RHIC, colhsional rather than 



o 

^^ ■ the radiative energy loss is the dominant mechanism for jet quenching. Ac- 



cordingly we calculate nuclear suppression factor for light hadrons by taking 



^~^' only the elastic energy loss in sharp contrast with the previous calculations 



where only the radiative loss are considered. 



X ■ PACS numbers: 12.38.Mh, 24.85.+p, 25.75.-q,13.87.Fh 



Jet quenching is one of the most promising tools to extract the initial parton density pro- 
duced in high energy heavy ion collisions. This is related to the final state energy loss of the 
leading partons [1-3] causing depopulation of hadrons at high transverse momentum (see [4] 
for experimental results). The suppressions of high pt hadrons and unbalanced back-to- 
back azimuthal correlations of the dijet events measured at Relativistic Heavy Ion Collider 
(RHIC) provide experimental evidence in support of the quenching. Based on the calcula- 
tions performed by several authors [5-7] the detailed theory of 'jet tomography' was devel- 
oped by Gyulassy et al. [3] considering only the energy loss due to induced bremsstrahlung 
radiation. The observed nuclear suppression of light hadrons (tt,-?]) in Au + Au collisions 
at ^/s = 62 — 200 AGeV at RHIC could be accounted for in these models. In all these 
analyses the collisional loss was ignored [8,9]. The non-photonic single electron spectrum 
from heavv meson decavs measured bv PHENIX Collaboration flOl nut this assumntion in 



question. No realistic parameter set can explain this data using the radiative energy loss 
based jet tomography model which either requires violation of bulk entropy bounds or non- 
perturbatively large Ug of the theory [12], or equivalently one requires excessive transport 
co-efficient geff = 14 GeVVfm [13]. 

The importance of coUisional loss in the context of RHIC was first discussed by the 
present authors [14,15]. It is shown in ref. [14] that there exists an energy range where 
collisional loss is as important as or even greater than its radiative counter part, hence 
cannot be neglected in any realistic model of jet quenching. Recently this is also noted in 
ref. [11,12,16,17]. It is similar to the passage of charged particles through material medium 
where the ionization loss is know to be the dominant mechanism at lower energies while at 
higher energies bremsstrahlung takes over. There exists a critical energy Ef. at which they 
contribute equally i.e. {dE/dx)rad = {dE/dx)cou aX E = Ec- For example, for an electron 
(proton) traversing copper target E^. ~ 25 MeV (1 GeV) [18]. Note that for heavier particle 
Ec is higher. This indicates that for the heavy quark collisional loss may be more important 
than the radiative loss at intermediate energies. In ref. [14] we have calculated Ec for light 
partons under RHIC conditions. 

In this light, we, in the present work would like to address if the omission of collisional 
loss at RHIC is justified or not. We argue that, whether the collisional or radiative loss 
is the main mechanism is a Pt dependent question. It also depends on the energies of the 
colliding system and expected to be different for RHIC and Large Hadron Collider (LHC). 
In contrast to the previous works, we, therefore, calculate nuclear suppression factor (Raa) 
for pions considering only the collision energy loss. At the end we shall show that there 
exists a pr window where this is reasonable assumption contrary to the commonly held view 
that collisional loss (for light partons) can be ignored altogether. 

The neutral pion production [19] (for charged hadrons see [20]) at RHIC in the pt 
window ~ 1 — 13 GeV, is found to be suppressed compared to the binary scaled p-p esti- 
mation [4]. This is attributed to the final state energy loss of the partons while passing 
through the nlasma, before frapmentinEr into hadrons f21— 231. The enerpv loss in the stan- 



dard perturbative calculations can be incorporated by modifying the fragmentation func- 
tion. This is accomplished by replacing fractional momentum z carried by the hadrons 
with z* = z/{l — Az) in the argument of the fragmentation function, D{z,Q'^), where 
Az = AE/E. This implementation assumes that all the partons suffer equal amount of 
energy loss which is questionable as argued in ref. [24,25]. We, therefore, take a different 
approach where the initial spectra is evolved dynamically by using Fokker Planck (FP) equa- 
tion. FP equation can be derived from Boltzmann equation if the collisions are dominated by 
the small angle scattering involving soft momentum exchange [15,26-31]. For an expanding 
plasma, FP equation takes the following form: 



+nirTlB±fiP,t)], (1) 
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where the second term on the left hand side arises due to expansion. Bjorken hydrodynamical 
model [32] has been used here for space time evolution. In Eq. (1) 1] denotes drag coefficient 
which is related to the energy loss or the 'stopping power' of the plasma, rj = {1/ E)dE/dx. 
B\\, Bj_ denote diffusion constants along parallel and perpendicular directions of the prop- 
agating parton representing rate of longitudinal and transverse broadening (variance) i.e. 
-B|| = d{{Ap\\)'^)/dt, B± = d{{Ap_[_Y) /dt. These transport coefficients can be calculated 
from the following expressions: 

dE ^ f d^kd^qdu 
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dx {2nf J 2k2k'2p2p' 

(-M)Lo/(|k|)[l + /(|k + q|)]a.. (2) 

^ u f d^kd^qduj ., , . , , 

(-M)Lo/(|k|)[l + /(|k + q|)]gl||. (3) 

In the above equations the small angle limit has been taken to write the arguments of the 
delta functions fl4.15l. 



The matrix elements include diagrams involving exchange of massless gluons which render 
7] and -B||,± infrared divergent. Such divergences can naturally be cured by using the hard 
thermal loop (HTL) [33] corrected propagator for the gluons as discussed below. We work 
in the coulomb gauge where the gluon propagator for the transverse and longitudinal modes 
are denoted by -Dqo = ^|| and Dij = {6ij — q^q^ /q'^)A± with [34]: 
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A±(go,g) ^ = ql-q'^ + 7:^1 
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The HTL modified matrix element in the limit of small angle scattering takes the fol- 
lowing form [14,15] for all the partonic processes having dominant small angle contributions 
hke qg ^' qg, qq —* qq etc. : 



\M\' = g'CRl6{EpEkY\A\\{qo,q) 
+ (vp X g).(vkX g)A^(go,g)|^ 



(6) 



where Cr is the appropriate color factor. With the screened interaction, the drag and 
diffusion constants can be calculated along the line of ref. [15]. 

Having known the drag and diffusion, we proceed to solve the FP equation. For this 
purpose we require the initial parton distributions which is parametrized as [35]: 



f{PT,Pz,t = ti) = 



dN 



Nn 
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where p^, a and Nq are parameters. Solving the FP equation with the boundary conditions, 
fiPi't) -^ for \p\ -^ oo, we are ready to evaluate the nuclear suppression factor, Raa 
defined as [36], 



Raa{P' 
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FIG. 1. Nuclear suppression factor for pion. Experimental data are taken from PHENIX col- 
laboration [19] for Au + Au collisions at ^/s = 200 GeV. Solid line indicates result from the 
present calculation with collisional energy loss of the partons propagating through the plasma 
before fragmenting into pions. 

where /(p',Tj) and /(p',Tc) denote the parton distributions at proper time Tj and Tc 
respectively. Here Tj is the initial time and Tc is the time when the system cools down to the 
transition temperature Tc (=190 MeV) [37]. The result for neutral pion is shown in Fig. 1 
which describes the PHENIX data [19] for Au + Au at ^/s = 200 GeV reasonably well. 

It should be noted here that the Raa{Pt) with collisional loss has a tendency to increase 
for higher pt, indicating less importance of collisional loss at this domain, where the radiative 
loss may become important. Therefore, a detailed calculation with both collisional and 
radiative loss may be useful to delineate the importance of individual mechanism. 
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FIG. 2. Excitation function of the nuclear modification factor for neutral pions in central A+A 
reactions at a fixed pr = 4 GeV where only the elastic energy loss is considered. Experimental 
data are taken from [36]. 

To stress our point further we also analyse the excitation function of the nuclear suppres- 
sion factor. Results are shown in Fig. 2. It is clear that Raa{pt) at pt = 4 GeV for various 
beam energies are found to be well described. The values of parton (quarks, anti-quarks and 
gluons) densities (^3+9+5) of the QCD medium which describe the data for various beam 
energies are shown in table I. 

To pin down the relative importance of 2 ^^ 2 2 -^ 3 processes, we determine the 
average energy of the parton which contribute to the measured Pt window of the hadrons. 
To this end, the average fractional momentum {{z)) of the fragmenting partons carried by 
the pion is calculated using relevant parton distribution and fragmentation functions. For 
the former, we use CTEQ [39] including shadowing via EKS98 parametrization [40], while 
for the fragmentation function KKP parametrization is used [41]. The average energy of the 
parton, (^E^^'^^"^) is obtained by using the relation < E^"''^^"^'' >= pJ^/ < z > for y^ = 0. The 
results are shown in Fig. 3. Our results are consistent with that of ref. [36] which quotes 
{z) = Phadron/Pparton - 0.5 - 0.7 for Phadron > 4 GcV at RHIC energies. 



TABLES 



17.3 4 0.97 

62.4 15 0.59 
130 27 0.39 
200 49 0.27 

TABLE I. The extracted initial parton densities and nuclear modification factors at various 

beam energies 
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FIG. 3. Average parton energy versus transverse momentum of pion for ^/s = 200 GeV/A. 



It might be recalled that at RHIC energies the nuclear modification factor Raa{pt) has 
been measured in the pion transverse momenta range pr ~ 1 — 13 GeV. Assuming that these 
pions are originated from the fragmenting partons, we ask the question, what is the average 
parton energy required to produce these pions? From Fig. 3, it is clear that the maximum 
average parton energy required is about 26 GeV here. 

Now the next question is, what is the dominant energy loss mechanism for partons with 
energy ~ 26 GeV or less? We might compare this value with the determined Ec^ given 
in refs. [14] (can also be read out from [16]) and note that at these energies collisional loss 
cannot be neglected. For lower beam energy, 62.4 (130) AGeV the value of maximum average 
parton energy required to produce a 13 GeV pion is 16 (22) GeV, where the collisional loss 
will definitely be more important. It is worthwhile to mention here that this estimation 
of Ec has some uncertainty as it depends on the length of the plasma, initial temperature, 
mean free path, dynamical screening mass etc. Those would affect both the mechanisms 
{i.e. radiative and collisional) of energy loss. Our chosen parameter set is consistent with 
that of ref. [7,42] used to study the radiative energy loss. 

In conclusion, our investigations clearly suggest that in the measured Pt range of light 



^Note that E^ is defined to be the energy below which elastic loss dominates [14]. 



hadrons at RHIC collisional, rather than the radiative, is the dominant mechanism of jet 
quenching. This is in sharp contrast to all the previous analyses. The determination of 
the critical energj{Ec), however, might change depending upon the detailed model of 'jet 
quenching'. Inclusion of three body elastic channels for heavy quark energy loss, which 
are considered in ref. [43], if applied for light flavours, might even increase E^, making our 
point stronger. E^ will also increase if there exist partonic bound state in the plasma due 
to ionization loss [44]. It should be mentioned that for the collisional energy loss we have 
not included finite size effect which however is shown to be small [16]. In light of these new 
findings the theory of jet tomography is expected to change considerably. 
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